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Abstract 

A systematic ab initio study of the second-order nonlinear optical properties of BN nanotubes 
within density functional theory in the local density approximation has been performed. Highly ac- 
curate full-potential projector augmented-wave method was used. Specifically, the second-harmonic 

(2) 

generation {x a h c ) an d linear electro-optical (r a b c ) coefficients of a large number of the single- 
walled zigzag, armchair and chiral BN nanotubes (BN-NT) as well as the double-walled zigzag 
(12,0)@(20,0) BN nanotube and the single-walled zigzag (12,0) BN-NT bundle have been calcu- 
lated. Importantly, unlike carbon nanotubes, both the zigzag and chiral BN-NTs are found to 

(2) 

exhibit large second-order nonlinear optical behavior with the X a b c an< ^ r abc coefficients being up 
to thirty times larger than that of bulk BN in both zinc-blende and wurtzite structures, indicating 
that BN-NTs are promising materials for nonlinear optical and opto-electric applications. Though 
the interwall interaction in the double-walled BN-NTs is found to reduce the second-order nonlin- 
ear optical coefficients significantly, the interwall interaction in the single-walled BN-NT bundle 
has essentially no effect on the nonlinear optical properties. The prominant features in the spectra 

(2) 

of X a b c (~ ui, to) of the BN-NTs are successfully correlated with the features in the linear optical 
dielectric function e{oj) in terms of single-photon and two-photon resonances. 
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I. INTRODUCTION 



Since their discovery in 1991 1 , carbon nanotubes (CNTs) have attracted considerable 
interest worldwide because of their unusual properties and great potentials for technological 
applications. For example, because of their one-dimensional character, metallic CNTs are 
quantum wires that may exhibit exotic Luttinger-liquid behavior rather than usual Fermi- 
liquid behavior in normal metal wires.— It was also predicted that nanotori formed from 
metallic CNTs may exhibit giant paramagnetic moments.— Furthermore, chiral CNTs are 
expected to exhibit a number of unusual optical properties such as optical activity, circular 
dichroism and second harmonic generation (see ^ and references therein) . 

Soon after the discovery of CNTs it became obvious that similar nanostructures could 
be formed by other elements and compounds which form layered structures bearing some 
resemblance to graphite. For example, hexagonal BN (/i-BN) was predicted on the basis of 
theoretical calculations^ to be capable of forming nanotubes, a prediction which was later 
confirmed experimentally by the synthesis of such nanotubes. 8 Both single- walled and multi- 
walled BN nanotubes (BN-NT) can now be readily synthesized? 9 Though CNTs continue to 
attract great interest, other nanotubes such as BN-NTs are interesting in their own right and 
may be able to offer different possibilities for technological applications that CNTs cannot 
provide. In particular, as far as the optical and opto-electronic applications of nanotubes 
are concerned, BN-NT could be superior to CNTs because BN-NTs are uniformly insulat- 
ing, independent of their chirality. Furthermore, BN-NTs tend to have a zigzag structure.— 
Though it is interesting that, depending on their chirality, CNTs can be metallic or semi- 
conducting or insulating^, it is still impossible to grow CNTs with a pre-specified chirality 
at present. Finally, recent experiments indicate that BN-NTs exhibit stronger resistance to 
oxidation at high temperatures than CNTs^i 

Therefore, the electronic, optical and other properties of both single-walled and mul- 
tiwalled BN-NTs are interesting and have been intensively studied theoretically in recent 
years (see, e.g., SiV 1 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ). In particular, Chen et a/.— calculated the transverse 
dielectric function of bundles of single-walled BN-NTs using a tight-binding model. Ng 
and Zhangi& calculated the optical absorption spectra of single-walled BN-NTs within a 
time-dependent localized-density-matrix approach based on a semiempirical Hamiltonian. 
Despite of these intensive theoretical studies, only few ab initio calculations of the optical 
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properties of BN-NTs have been reported— because of the heavy demand of the computing 
resources. Semiempirical tight-binding model is known to describe well only the electronic 
excitations near the band gap of the large radius BN-NTs. Systematic ab initio calculations 
of the optical properties are thus needed in order to quantitatively interpret the optical 
experiments and to predict the BN-NTs with desired optical properties. Therefore, we have 
recently carried out a series of ab initio calculations^ in order to analysize the band struc- 
ture and linear optical features of all the three types of the BN-NTs and their possible 
dependence on diameter and chirality. In this work, we investigate the second-order op- 
tical susceptibility and also linear electro-optical coefficient of the BN-NTs. The primary 
objective is to find out the features and magnitude of the second-harmonic generation and 
linear electro-optical coefficients of the BN-NTs in order to see whether they have any po- 
tential applications in nonlinear optical and electro-optical devices such as second-harmonic 
generation, sum-frequency generation and electrical optical switch. The second objective 
is to identify characteristic differences in nonlinear optical properties between BN-NTs and 
CNTs. 

The rest of this paper is organized as follows. In Sec. II, the theoretical approach and 
computational details are briefly described. In Sec. Ill, the calculated second-order nonlinear 
optical susceptibility and linear electro-optical coefficients of /i-BN, single hexagonal BN 
sheet, and BN nanotubes are presented and analysized. Finally, in Sec. IV, a summary is 
given. 

II. THEORY AND COMPUTATIONAL METHOD 

BN-NTs can be considered as a layer of graphitic BN sheet rolled up into a cylinder, 
and the structure of a BN-NT is completely specified by the chiral vector which is given in 
term of a pair of integers (n,m).^l As for CNTs, BN-NTs are classified into three types, 
namely, armchair (n,n) nanotubes, zigzag (n,0) nanotubes, and chiral (n,m) nanotubes with 
n ^ m.^2i We consider a large number of representative BN-NTs with a range of diameters 
from all three types, namely, the zigzag [(n,0), n — 5, 6, 8, 9, 12, 13, 15, 16, 17, 20, 21, 24, 
25, 27], armchair [(n,n), n = 4, 5, 6, 8, 12, 15], and chiral (4,2), (6,4), (8,4), (10,5) BN-NTs. 
The double-walled zigzag (12,0)@(20,0) BN nanotube and the single-walled zigzag (12,0) 
BN-NT bundle have also been investigated to see the effects of interwall interaction. Our 
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ab initio calculations for the BN-NTs were performed using highly accurate full-potential 
projector augmented-wave (PAW) method 2 ^, as implemented in the VASP package 2 ^. They 
are based on density functional theory (DFT) with the local density approximation (LDA). 
A supercell geometry was adopted so that the nanotubes are aligned in a square array with 
the closest distance between adjacent nanotubes being at least 6 A. A large plane- wave 
cut-off of 450 eV was used throughout. 

Firstly, the ideal nanotubes were constructed by rolling- up a hexagonal BN sheet. Their 
atomic positions and lattice constants were then fully relaxed by a conjugate gradient tech- 
nique. Theoretical equilibrium nanotube structures were obtained when the forces acting on 
all the atoms and the uniaxial stress were less than 0.03 eV/A and 2.0 kBar, respectively. 
The theoretical equilibrium lattice constants T and curvature energies E c (total energy rel- 
ative to that of single BN sheet) as well as the computational details have been reported 
before 2 ^ 

The self-consistent electronic band structure calculations were then carried out for the 
theoretically determined BN-NT structures. In this work, the non-linear optical properties 
were calculated based on the independent-particle approximation, i.e., the excitonic effects 
and the local-field corrections were neglected. As reported before 2 ^, the dielectric function 
of /i-BN calculated within the single-electron picture are in reasonably good agreement with 
experiments. Therefore, it might be expected that the independent-particle approximation 
could work rather well for the non-linear optical properties of the BN-NTs too. Following 
previous nonlinear optical calculations 4 * 24 , the imaginary part of the second-order optical 
susceptibility due to direction interband transitons is given by 2 ^ 

where the contribution due to the so-called virtual-electron (VE) process is 
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and that due to the virtual-hole (VH) process 



5(e i; - w) + — 3 <Heji -2a;)}. (3) 



Here e^j = e^j — e^i and (p^p^) = ^(P^iPu +PuP C ji)- The dipole transition matrix elements 
P\a =< kj|j3 a |ki > were obtained from the self-consistent band structures within the PAW 
formalism—. The real part of the second-order optical susceptibility is then obtained from 
XabJ by a Kramer-Kronig transformation 



X [ >{-2u,u,u) = -P / du; - 



7T /,, Co>' 2 — Co> 2 



(4) 



The linear electro-optic coefficient r abc {u) is connected to the second-order optical sus- 
ceptibility Xabc(~ u i 0) through the relation 2 ^ 

xi? c (-^, 0) = ~n 2 a {uj)n 2 b {uj)r abc {uj) (5) 

where n(u) is the refraction index in the a-direction. In the zero frequency limit, 



Therefore 



lin J,xi 6 j c (-2a;,a;,a;) = lim lx^-w, w, 0). (6) 



Furthermore, for the photon energy ftu; well below the band gap, the linear electro-optic 
coefficient r abc (u) ~ r abc (0) because Xa& c ( — 2u;, a>, cu) and n(w) are nearly constant in this low 
frequency region, as shown in the next Sec. and in Ref. — . 

In the present calculations, the 5-function in Equs. 2-3 is approximated by a Gaussian 
function with T = 0.2 eV. The same fc-point grid as in the DOS calculation is used. Further- 
more, to ensure that e' calculated via Kramer-Kronig transformation (Equ. 4) is reliable, at 
least ten energy bands per atom are included in the present optical calculations. The unit 
cell volume Q in Equs. 2-3 is not well defined for nanotubes. Therefore, like the previous 
calculations^ 2 ^, we used the effective unit cell volume of the nanotubes rather than the 
volume of the supercells which is arbitrary. The effective unit cell of a nanotube is given by 
Q = tt[(D/2 + d/2) 2 - (D/2 - d/2) 2 ]T = ixDdT where d is the thickness of the nanotube 
cylinder which is set to the interlayer distance of /i-BN (3.28 A 2 ^). D and T are the diameter 
and length of translational vector of the nanotube 20 , respectively. 



III. RESULTS AND DISCUSSION 



A. Hexagonal BN and single BN sheet 

For comparison with the BN-NTs, we first investigated the second-order nonlinear optical 
susceptibility of both /i-BN and an isolated honeycomb BN sheet. The isolated BN sheet is 
simulated by a slab-supercell approach with an inter-sheet distance of 6.5 A. The theoret- 
ically determined lattice constants (a = 2.486 A and c = 6.562 A for /i-BN and a = 2.485 
A for the BN sheet were used. Note that the theoretical lattice constants of /i-BN agree 
rather well (within 1.5 %) with the experimental values (a = 2.50 A and c = 6.65 A)2£. 
Interestingly, we find numerically that although /i-BN has zero second-order nonlinear opti- 

(2) (2) (2) 

cal susceptibility, the Xaabi Xbaa an d Xm f° r the isolated BN sheet are nonzero. Here a and 

(2) (2) 

b denote the two Cartesian coordinates within the BN layer. Moreoever, xla a = Xaab ano - 
Xbbb = ~Xaab- This is consistent with the symmetry consideration, demonstrating that our 
numerical method and calculations are qualitatively correct. Bulk /i-BN should have zero 
second-order nonlinear optical susceptibility due to its spatial inversion symmetry. On the 
other hand, the isolated BN sheet does not have the spatial inversion symmetry (Dg h ) and 
its symmetry class is D 3 h (P6ra2). Therefore, the isolated BN sheet has nonzero Xaabi Xbaa 
and Xbbb fr° m the symmetry consideration. 29 

In Fig. 1, the calculated real and imaginary parts as well as the absolute value of the 
imaginary part of Xmb(~ < ^ UJ -> u i u ) °f the single BN sheet are shown. In the calculations, a 
/c-point grid of 100 x 100 x 2 for the honeycomb BN sheet is used. It is clear from Fig. 1 
that the second harmonic generation (SHG) coefficient Xaab(~ ^ ^0 ^ s significant in the 
entire range of the optical photon energy (hw). Furthermore, for the photon energy smaller 
than 2.2 eV, the x a ab * s P ure ly dispersive (i.e., real and lossless) (Fig. la), suggesting that 
the BN sheet has potential application in nonlinear optical devices. Note that since the BN 
sheet has a theoretical band gap (E g ) of ~4.5 eV, the absorptive (imaginary) part of the 
X^alb becomes nonzero only above the half of the band gap (i.e., ~2.2 eV). The real part of 
the Xaab remains nearly constant at low photon energies up to 1.0 eV, then increases steadily 
in magnitude as the photon energy increases, and finally peaks at the absorption edge of 
~2.2 eV (Fig. la). In the energy range from 2.5 to 5.5 eV, the real part of the Xaab becomes 
positive and forms a broad double peak structure. Beyond 5.5 eV, it becomes negative again 
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FIG. 1: (a) Real and imaginary parts as well as (b) the absolute value of the imaginary part of 
X a ^ b (—2u> , u> , uj) of the isolated BN sheet. In (c), s"(u>) and e"(u>/2) (imaginary part of the dielectric 
function) from Ref. — are plotted. 

and its magnitude gradually deminishes as the photon energy further increases (Fig. la). 

(2) 

The absorptive part of the Xaab 1S nonzero only above ~2.2 eV, and looks like a Lorentzian 
oscillation between 2.2 and 6.0 eV with one sharp negative peak at ~ 2.5 eV and one broad 
positive peak around 5.5 eV (Fig. la). It is clear from Equs. 2-3 that the calculated Xaab 
spectra can have pronounced features due to both single- and double-frequency resonant 
terms. To analyze the features in the calculated spectra, it is helpful to compare the 
absolute value of x"^ (Fig lb) with the absorptive part of the corresponding dielectric 
function e"M Therefore, the calculated e" from our previous publication^ are shown in 
Fig. lc as a function of both uj/2 and u. Clearly, the first sharp peak at ~2.5 eV is due 
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TABLE I: Calculated static refraction index n, second-order optical susceptibility (0) and linear 
electro-optical coefficient r abc of the isolated BN sheet. 



n a (n c ) 


Xm, Xaab (P m / V ) 


nbb, r aab (pm/V) 


2.19 (1.63) 


81.3, -81.3 


-3.56, 3.56 



to two-photon resonances [cf. e"(u>/2)] whilst in contrast, the second broad peak around 
5.5 eV comes from the single-photon resonances [cf. e^(a>)]. Nevertheless, both single- and 
double-photon resonances involve only interband 7r — > it* and a — > a* optical transitions for 
the electric field vector E polarized parallel to the BN layer (E || a)^. 

In Table I, the calculated zero frequency linear electro-optic coefficient r(0) as well as 
the corresponding second-order nonlinear optical susceptibility x^ (0,0,0) and the static 
refraction index n(0) are listed. The refraction index n(0)(= a/^(0)) is derived from the 
calculated static dielectric constant e(0) which has been reported in our recent publication^. 
Note that the r(0) and x (2) (0, 0, 0) are listed in the SI pm/V unit, and 1 pm/V = 4. 1888 x 1(T 8 
esu. Significantly, the static second-order optical susceptibility for the isolated BN sheet is 
nearly thirty times larger than that of BN in both the zinc-blende and wurtzite structures^. 

B. Second-order optical susceptibility 

We have explicitly calculated the second-order optical susceptibility for the zigzag (5,0), 
(6,0), (8,0), (9,0), (12,0), (13,0), (15,0), (16,0), (17,0), (20,0), (21,0), (24,0), (25,0), (27,0), 
armchair (3,3), (4,4), (5,5), (6,6), (8,8), (12,12), (15,15), and chiral (4,2), (6,4), (8,4), (10,5) 
BN-NTs. In the calculations, a uniform grid (1 x lxm) along the tube axis (z-axis) was used. 
The number m is 200, 120, 150, 60, 100, 100, 80, 120, 80, 80, 60, 80, 60, 60 for the sigzag 
BN-NTs, respectively, and 60, 50, 40, 60 for the chiral BN-NTs, respectively, In the case 
of CNTs, only the chiral nanotubes would exhibit second-order nonlinear optical behavior 
with two nonvanishing components of xyz and xzy of x*^-' 4 Here z refers to the coordinate 
along the tube axis whilst x and y denote the two coordinates that are perpendicular to 
the tube axis. As for CNTs, the armchair BN-NTs are found not to have any nonzero 
components of x^ ■ I* 1 contrast, both the zigzag and chiral BN-NTs are found to show 
second-order nonlinear optical behavior. Specifically, all the zigzag BN-NTs except (5,0), 
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(i) (6,0) BN nanotube (ii) (12,0) BN nanotube (iii) (20,0) BN nanotube (iv) (25,0) BN nanotube 
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FIG. 2: (Color online) Real and imaginary parts of x^( — 2tv, u, ui) of the zigzag (6,0), (12,0), 
(20,0) and (25,0) BN nanotubes. 

(9,0) and (27,0), have six nonvanishing components of the second-order optical susceptibility, 
namely, xzx, xxz, yyz, zxx, zyy, zzz. Nevertheless, these components are not completely 



independent of each other. In particular, Xxiz = Xvvz = Xx%, and xfyy = X^xx- This finding 



,(2) 



.(2) 



is especially important for the application of BN-NTs in nonlinear optical devices because 
most BN-NTs tend to have a zigzag structure^ as mentioned before. The chiral BN-NTs have 
eight nonvanishing components of the second-order optical susceptibility with two additional 

(2) (2) 

nonzero components being xyz and yzx. Note that Xyzx — — Xxyz- These numerical findings 
are consistent with the consideration of the symmetry of the BN nanotubes. The point 
symmetry groups of the BN-NTs^ are Ci nv for zigzag (n,0) nanotubes, Cmh fo r armchair 
(n,n) nanotubes, and Cn for chiral (n,m) nanotubes where N = 2(n 2 + m 2 + nm)/dR with 
cIr being the greatest common divisor of 2n + m and 2m + n. Therefore, these symmetries 
would dictate^ that all the components vanish for armchair (n, n) nanotubes, and that 
nonvanishing components for zigzag (n,0) nanotubes are xzx = yzy, xxz = yyz, zxx = zyy, 
zzz, as well as that nonvanishing components of chiral (n, m) nanotubes include all that 
for zigzag (n,0) nanotubes plus xyz = -yzx. We don't know at the moment why the (5,0), 
(9,0) and (27,0) BN-NTs are found numerically to have no nonzero components though they 
don't have to from their symmetry point of view. 

In Fig. 2, the calculated real and imaginary parts of the second-order optical suscep- 
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tibility x (2) (- 2 ^, uj,uj) for the four selected zigzag nanotubes [(6,0), (12,0), (20,0), and 
(25,0)] are shown. As for the single BN sheet, the second harmonic generation coefficients 
t^( 2 )(— 2u, u, uS) for the zigzag BN-NTs are significant in the entire range of the optical pho- 
ton energy (fuv). Indeed, they are generally more than ten times larger than that of bulk 
BN in the zinc-blende and wurtzite structures.— Moreover, for the photon energy smaller 
than 2.0 eV, the x^ is purely dispersive (i.e., real and lossless) (Fig. 2). Interestingly, the 

(2) 

magnitude of Xzzz is the largest and in general about three times larger than that of the 

(2) 

other nonvanishing components (Fig. 2). For Xzzz, the electric field of both the incoming 
and outgoing photons is polarized parallel to the tube axis and thus the electric dipolar- 
ization effect would be essentially zero^S. This may be particularly important for nonlinear 

(2) (2) 

optical applications. We also note that the shape of the spectra of Xzzz an d Xzxx for all 
the zigzag BN-NTs look very similar, except the difference in sign. For the zigzag BN-NTs 

(2) 

with a larger diameter such as (12,0), (20,0) and (25,0), the spectrum of XxJx is also similar 
to that of x^zxx (Fig. 2). The magnitude of all the components decreases somewhat as the 
diameter of the tubes increases [e.g., from (6,0) to (12,0)], but however, becomes stable as 
the diameter further increases [see, e.g., (20,0) and (25,0) in Fig. 2]. Furthermore, both the 

(2) 

magnitude and shape of the Xzzz spectrum for the zigzag BN-NTs with a larger diameter, 
e.g., (25,0), approach to that of the single BN sheet (Fig. 1 and Fig. 2), as they should. 

As for the single BN sheet, in order to understand the features in the calculated x 
spectra of the zigzag BN-NTs, the absolute values of the imaginary part x"^ °f a H the 
nonzero components of the (12,0) BN-NT are plotted and compared with the absorptive 
part of the dielectric function e" from our previous publication 20 in Fig. 3. Strikingly, the 
first prominent peak between 2.0 and 4.0 eV in the Xzzz spectrum is almost identical to the 
first peak in the e" z {uj/2) (see Fig. 3d-e), indicating that it is due to two-photon resonances. 
In contrast, the second peak between 4.5 and 5.5 eV in the Xzzz spectrum is very similar 
to the first peak in the e"(u)), suggesting that it is caused by the single-photon resonances. 
Nevertheless, both these single- and double-photon resonances involve only interband tt — > it* 
and u — ► <7* optical transitions for the electric field vector E polarized parallel to the tube 
axis (E || z) 20 . Fig. 3 also shows that the double-peak structure between 2.0 and 4.0 eV 
in both the x'xzx and x'zix spectra is mainly due to the two-photon resonances with E _L z 
[cf. e"(u/2)] (see Fig. 3a-c), while, in contrast, the second feature in the photon energies 
above 4.5 eV perhaps comes predominantly from the the single-photon resonances for E _L z 
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FIG. 3: Absolute value of the imaginary part of x^ 2 \~ 2<^,<^,<^) (a, c, d) as well as s'^(u) and 
e"(to/2) (imaginary part of the dielectric function) (b, e) from Ref. — of the zigzag (12,0) BN 
nanotube. 

[cf. e^(u;)]. This conclusion is further supported by the fact that the magnitude of s"(uj) 
(e"(uj/2)) is only about half of that of s"(oo) (e"(ui/2)), and concurrently the magnitude of 
Xxzx and x'zxx is only about half of that of x'ziz too. 

In Fig. 4, the calculated real and imaginary parts of the second-order optical susceptibility 
X^ 2 \—2u,u,uj) for all the four chiral nanotubes [(4,2), (6,2), (8,4) and (10,5)] are displayed. 
In general, the spectra of each component of the second-order optical susceptibility for all the 
chiral BN-NTs are similar, except that the magnitude decreases somewhat as the diameter 
of the tubes increases. The decrease of the magnitude with the tube diameter is particularly 
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(i) (4,2) BN nanotube (ii) (6,2) BN nanotube (iii) (8,4) BN nanotube (iv) (10,5) BN nanotube 
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FIG. 4: (Color online) Real and imaginary parts of x^ 2 \~ 2w,cj, oS) of the chiral (4,2), (6,2), (8,4) 
and (10,5) BN nanotubes. 

apparent for the xyz and zzz components. For all the chiral BN-NTs, remarkably, the zzz 
component is nearly ten times larger than all the other nonvanishing components. Another 
common feature is that both the real and imaginary parts of x^ 2 \~ 2a;, to, 10) show a rather 
oscillatory behavior, particularly for the xyz and xzx components (Fig. 4). The amplitude 
of these oscillatory real and imaginary parts is rather large in the photon energies of 2.0~ 
6.0 eV for the chiral BN-NTs with a small diameter. It should also be noted that the shape 
and magnitude of the xzx and zxx components look very much alike (Fig. 4b-c). 

Again, in order to understand the structures in the calculated x^ spectra of the chiral 
BN-NTs, the absolute values of the imaginary part x! °f an the nonzero components of 
the (6,2) BN-NT are plotted and compared with the absorptive part of the corresponding 
dielectric function e" from our previous publication 2 ^ in Fig. 5. Remarkably, the first 
prominent peak between 2.0 and 4.0 eV in the x'zzz spectrum is almost identical to the first 
peak in the e" z (ui/2) (see Fig. 5d-e), indicating that it is due to two-photon resonances. On 
the other hand, the second structure between 4.5 and 5.5 eV in the x'zzz spectrum may 
be correlated with the first peak in the e"(u), suggesting that it is caused by the single- 
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photon resonances. As for the zigzag BN-NTs, both these single-photon and double-photon 
resonances involve only interband 7r — > tt* and a — > a* optical transitions for the electric 
field vector E polarized parallel to the tube axis (E \\ z)^. Fig. 5 further suggests that the 
feature between 2.0 and 4.0 eV in the spectra of both the x'xzx and x'zxx as well as x'xyz may 
be attribued to the two-photon resonances with E _L z [cf. e"(uj/2)} (see Fig. 5a-c), while, 
in contrast, the second feature in the photon energies above 4.5 eV is mainly due to the 
single-photon resonances for E _L z [cf. e"(uj)}. One exception is that the narrow peak at 
2.5 eV in the x'zxx spectrum (Fig. 5c) which may come predominantly from the two-photon 
resonances with E || z [cf. e" z {uj/2) in Fig. 5e]. 



13 



TABLE II: Calculated static refraction index n, second-order optical susceptibility x^ an d linear 
electro-optical coefficient of the zigzag and chiral BN nanotubes. 
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1.95 (2.20) 
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-1.10, -1.12, 1.65 


(12,0) bundle 


1.92 (2.19) 


37.0, 38.3, -89.0 


-4.19, -4.33, 7.74 


(4,2) 


1.92 (2.14) 


11.5, 11.5, -61.8 


-1.36, -1.36, 5.89 


(6,2) 


1.89 (2.15) 


15.8, 15.8, -64.9 


-1.91, -1.91, 6.07 


(8,4) 


1.89 (2.15) 


16.0, 16.0, -40.6 


-1.94, -1.94, 3.80 


(10,5) 


1.89 (2.15) 


17.5, 17.5, -41.3 


-2.12, -2.12, 3.87 



C. Linear electro-optical coefficient 

In Table II, the calculated zero frequency linear electro-optic coefficient r(0) as well as 
the corresponding second-order nonlinear optical susceptibility x^(0, 0, 0) are listed. The 
r(0) is calculated from the corresponding (0,0,0) by using Equ. 7. The refraction 
index n(0)(= a/^(0)) is derived from the calculated static dielectric constant 5(0) which 
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has been reported in our recent publication^ .. Table II shows that apart from the (5,0), 
(9,0) and (27,0) BN-NTs which have no nonvanishing r a b c (0), all the other zigzag BN-NTs 
have very similar linear electro-optical coefficients, as for the static dielectric constant and 
polarizability 20 . Nevertheless, r zzz (0) decreases slightly as the diameter increases, while 
r xzx(fy and r zxx (0) increase slightly as the diameter goes up. We note that Xxzx and xixx 
are virtually identical, while Xzzz is about two times larger than Xxzx and xiaL- 

(2) (2) 

The chiral BN-NTs have two additional nonvanishing components Xxyz and Xyzx- Never- 

(2) (2) 

theless, the calculated static values of Xxyz and Xyzx are zero, satisfying the requirement by 
the so-called Kleinman symmetry^ which demands that x*yz(0) = Xyz\{p>)- Consequently, 
the corresponding static linear electro-optical coefficients r xyz (0) and r yzx (0) are zero too. 
Therefore, Xyzx, Xxyz, r xyz (0) and r yzx (0) for the chiral BN-NTs are not listed in Table II. 
As for the nonvanishing static components, r xzx (0) and r zxx (0) for the (6,2), (8,4) and (10,5) 
BN-NTs are rather close, while r zzz (0) for the (6,2) BN-NT is nearly 1.6 times larger than 
that for the (8,4) and (10,5) BN-NTs. Though r zzz (0) for the (4,2) BN-NT is close to that 
of the (6,2) BN-NT, r xzx (0) and r zxx (0) for the (4,2) is somewhat smaller than that of the 
other armchair BN-NTs. 

Group III nitrides are promising materials for nonlinear optical and opto-electric appli- 
cations because of their unique physical properties such as a wide band gap. However, the 
static X^(0) an d hence r(0) of bulk BN structures are the smallest among the group III 
nitrides (see, e.g.,—). It is therefore remarkable that BN in the single- walled zigzag BN-NTs 
and their bundles has greatly enhanced static second-order nonlinear optical susceptibility 
and linear electro-optical coefficient. This enhancement can be as large as thirty folds (see 
Table II and Ref. — ), and it could make the BN nanotube structures have the largest x^ 
and r among the nitrides. In fact, the static x (0) °f the zigzag BN-NTs is up to five times 
larger than that of GaN which has the largest x (0) among the bulk group III nitrides.— 

D. Effects of interwall interactions 

We have so far focused only on the nonlinear optical properties of single-walled BN-NTs. 
However, BN-NTs are usually multiwalled or in the format of bundles. Therefore, to study 
the possible effects of interwall interactions on the nonlinear optical properties of BN-NTs, 
we have also considered a double-walled BN-NT, namely, the zigzag (12,0)@(20,0), and a 
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FIG. 6: (Color online) Real and imaginary parts of X ( — 2w, uj, uj) of the double-walled 
(12,0)@(20,0) BN nanotube. 

single-walled zigzag (12,0) BN-NT bundle. The double-walled nanotube is chosen because 
BN-NTs typically have a zigzag structure, and also the interwall distance of about 3.2 A 
between the (12,0) and (20,0) nanotubes is close to the interlayer distance in /i-BN. The 
BN-NT bundle was modelled by a two-dimensional hexagonal array with the initial interwall 
distance between adjacent nanotubes being about 3.2 A. The atomic positions and lattice 
constants were then fully relaxed by a conjugate gradient technique. Theoretical equilibrium 
nanotube structures were obtained when the forces acting on all the atoms and the uniaxial 
stress were less than 0.03 eV/A and 2.0 kBar, respectively. The theoretical determined 
atomic structures were used in the electronic structure and optical property calculations. 
The calculated density of states and dielectric function of the (12,0)@(20,0) double-walled 
BN-NT have already been reported in Ref. — . 

The calculated second-order optical susceptibility of the double-walled (12,0)@(20,0) BN- 
NT is shown in Fig. 6. A uniform grid of 1 x 1 x 60 in the Brillouin zone was used. Figs. 2 
and 6 show that broadly speaking, each component of the second-order optical susceptibility 
of the (12,0)@(20,0) BN-NT is rather similar to the corresponding component of the (12,0) 
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and (20,0) BN-NTs. However, the magnitude of the second-order optical susceptibility of 
the double-walled BN-NT is significantly reduced in comparison with that of the (12,0) and 
(20,0) BN-NTs. In particular, the magnitude of the zzz component of the (12,0)@(20,0) 
BN-NT is about three and two times smaller than that of the (12,0) and (20,0) BN-NTs, 
respectively (Figs. 2 and 6). In contrast, the density of states and dielectric function of the 
(12,0)@(20,0) as well as (12,0) and (20,0) BN-NTs are very similar,^ indicating that the 
effects of interwall interaction on the electronic structure and linear optical properties would 
be small. In particular, the dielectric function of the (12,0)@(20,0) is nearly the same as that 
of the (20,0). The reduction in the magnitude of x ma Y be explained by the fact that bulk 
/i-BN does not show second-order nonlinear optical behavior and that a multi-walled BN-NT 
with a very large diameter is essentially equivalent to /i-BN. Another discernable difference 
is that the second-order optical susceptibility of the (12,0)@(20,0) BN-NT is considerably 
more oscillatory than that of the (12,0) and (20,0) BN-NTs. 

The static SHG and linear electro-optic coefficients are compared with the single-walled 
BN-NTs in Table II. Although the static refraction index of the (12,0)@(20,0) BN-NT is 
almost identical to that of the single-walled BN-NTs, the static SHG coefficient is more 
than three times smaller than that of the single-walled zigzag nanotubes. Consequently, the 
linear electro-optical coefficient of the double-walled nanotube is several times smaller than 
that of the single- walled zigzag nanotubes (Table II). 

The calculated second-order optical susceptibility of the (12,0) BN-NT bundle is shown 
in Fig. 7. A uniform grid of 6 x 6 x 60 in the Brillouin zone was used. The effective unit cell 
volume is used in order to compare quantitatively the optical properties of the bundle with 
that of the isolated (12,0) BN nanotube. Nevertheless, the use of the solid unit cell volume 
of the bundle would only uniformly reduce the amplitude of the calculated second-order 
optical susceptibility by a factor of 1.43 but would not change its shape. Figs. 2 and 7 
show that each component of the second-order optical susceptibility of the (12,0) BN-NT 
bundle is similar to the corresponding component of the isolated single- walled (12,0) BN-NT 
in both the shape and magnitude. In particular, the zzz component of the (12,0) BN-NT 
bundle is nearly identical to that of the isolated (12,0) BN-NT (Figs. 2 and 7). This clearly 
shows that the interwall interaction between the BN-NTs in the bundle has essentially no 
effect on the second-order nonlinear optical properties, in strong contrast to the case of the 
multi-walled BN-NTs. We also find that the dielectric function of the (12,0) BN-NT bundle 
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FIG. 7: (Color online) Real and imaginary parts of x^ 2 \~ 2uj,uj,uj) of the single-walled (12,0) 
BN-NT bundle. 

(not shown here) is nearly the same as that of the isolated single-walled (20,0) BN-NT. 
Nevertheless, there are minor differences in the second-order optical susceptibility between 
the isolated (12,0) BN-NT and the (12,0) BN-NT bundle. For example, the second-order 
optical susceptibility of the isolated (12,0) BN-NT is somewhat more oscillatory than that 
of the (12,0) BN-NT bundle. This minor difference might be due to the fact that more 
Appoints for the (12,0) BN-NT bundle were used than for the single (12,0) BN-NT. 

The static SHG and linear electro-optic coefficients are compared with the single-walled 
BN-NTs in Table II. Unlike the case of the double walled BN-NT, the static SHG linear 
and electro-optic coefficients of the single-walled (12,0) BN-NT bundle are rather similar 
to that of the single- walled BN-NTs (Table II). This clearly shows that the interwall in- 
teraction hardly has any effect on the nonlinear optical properties, in strong contrast to 
the double walled BN-NT where the interwall interaction significantly reduce the SHG and 
linear electro-optic coefficients. Indeed, the SHG and linear electro-optic coefficients of the 
(12,0) BN-NT bundle are slightly enhanced, as compared with that of the single (12,0) 
BN-NT (Table II). This is rather significant because in the nonlinear optical applications, 
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one would usually need the BN-NTs in the format of either multiwalled BN-NTs or BN-NT 
bundles. The present work therefore suggests that the single-walled BN-NT bundles would 
be prefered for nonlinear optical and electro-optical applications. 

IV. SUMMARY 

We have carried out a systematic ab initio study of the second-order nonlinear optical 
properties of the BN-NTs within density functional theory in the local density approxi- 
mation. We used the highly accurate full-potential PAW method. The underlying atomic 
structure of the BN nanotubes was determined theoretically. Specifically, the properties of 
the single-walled zigzag [(5,0), (6,0), (8,0), (9,0), (12,0), (13,0), (15,0), (16,0), (17,0), (20,0), 
(24,0), (25,0), (27,0)], armchair [(3,3), (4,4), (5,5), (6,6), (8,8), (12,12), (15,15)], and chiral 
[(4,2), (6,2), (8,4), (10,5)] nanotubes as well as the double-walled (12,0)@(20,0) nanotube 
and the single-walled (12,0) nanotube bundle have been calculated. For comparison, the 
second-order nonlinear optical properties of /i-BN and the single BN sheet have also been 
calculated. Interestingly, we find that though /i-BN has zero second-order nonlinear optical 

(2) (2) (2) 

susceptibility, the Xaabi Xbda anc ^ Xm f° r ^ ne isolated BN sheet are large and generally several 
ten times larger than that of bulk BN in both the zinc-blende and wurtzite structures. We 
also find that, unlike carbon nanotubes, both the chiral and zigzag BN-NTs have pronounced 
second-harmonic generation and linear electro-optical coefficients which are comparable to 
that of the single BN sheet. The prominant structures in the spectra of 2u, u, uo) 

of the BN-NTs have been successfully correlated with the features in the corresponding 
linear optical dielectric function s{uj) in terms of single-photon and double-photon reso- 
nances. Though the interwall interaction in the double-walled BN-NTs is found to reduce 
the second-order nonlinear optical coefficients significantly, the interwall interaction in the 
single-walled BN-NT bundle has essentially no effect on the nonlinear optical properties, 
suggesting that the single-walled BN-NTs and single-walled BN-NT bundles are promising 
nonlinear optical materials for applications in, e.g., second-harmonic generation, sum fre- 
quency generation and electro-optical switches. We hope that this work would stimulate 
experimental investigations into the second-order nonlinear optical properties of BN-NTs. 
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